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Abstract 

We present a detailed calculation of next-to-leading order QCD corrections to the Wt associated 
production using the one cutoff phase space slicing method. Such QCD corrections have been 
calculated independently by two groups already, however, a number of differences were found. It 
is desirable to have a third party calculation to make a crossing check. In this note, we present 
our complete results of the virtual corrections which are not shown in the literature so far. The 
numerical comparison will be presented in the forthcoming paper. As a demonstration of the 
one cutoff phase space slicing method, we also show in details how to organize the color ordered 
amplitudes and how to slice the soft and collinear phase space regions. 
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I. INTRODUCTION 
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There are three separate single top quark production processes of interest at the hadron 
collider, which may be characterized by the virtuality of the W boson (of four momentum 
q) in the processes. The s-channel process qq' — > W* — > tb via a virtual s-channel W boson 
involves a timelike W boson, q 2 > (m t + m&) 2 , the t-channel process qb — > q't (including 
q'b — > qt, also referred as IV-gluon fusion) involves a spacelike W boson, q 2 < 0, and the 



tW associated production process bg — > tW involves an on-shell W boson, q 2 



ni 



w 



Therefore, these three single top quark production mechanisms probe the charged-current 
interaction in different q 2 regions and are thus complementary to each other. To improve 
the theory prediction on the single top production rate, the next-to-leading order (NLO) 



corrections, at the order of a s , for these three channels has been carried out in Refs. 
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The NLO QCD corrections to the tW associated production have been calculated inde- 
pendently by two groups 
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35| , but a number of differences were found. It is desirable to 



have a third party calculation to make a crossing check. Furthermore, there is no analytic 
result available in the literature so far. In this note we present a detailed calculation of the 
NLO QCD corrections to the tW associated production using one cutoff phase space slicing 
method. The rest of this paper is organized as follows. In Sec. II, we briefly review the 
method of our calculation. In Sec. Ill, we present the Born level helicity amplitudes of the 
tW production. In Sec. IV we present the NLO virtual corrections. In Sec. V we present 
the calculation of the soft and collinear singularities of the real emission correction using 
the one cutoff phase space slicing method. 



II. SHORT REVIEW OF THE PHASE SLICING METHOD 

The construction of a flexible event generator requires the generation of partonic final 
states with a minimal amount of implicit phase space integration. At the leading order this 
is trivial, but in the calculation of NLO QCD corrections, one generally encounters both 
ultraviolet (UV) and infrared (IR) (soft and collinear) divergences. The former divergences 
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can be removed by proper renormalization of couplings and wave functions. In order to 
handle the latter divergences, one has to consider both virtual and real corrections and 
carefully handle the cancellation of divergences between the soft and collinear contributions 
and the virtual corrections. The soft divergences will cancel according to the Kinoshita- 
Lee-Nauenberg (KLN) theorem 3y, [37|, but some collinear divergences remain uncancelled. 
In the case of considering the initial state partons, one needs to absorb additional collinear 
divergences to define the NLO parton distribution function of the initial state partons. After 
that, all the infrared-safe observables will be free of any singularities. To calculate the inclu- 
sive production rate, one can use dimensional regularization to regularize divergences and 
adopt the modified minimal subtraction (MS) factorization scheme to obtain the total rate. 
However, owing to the complicated phase space for multi-parton configurations, analytic 
calculations are in practice impossible for all but the simplest quantities. During the last 
few years, effective numerical computational techniques have been developed to calculate 
the fully differential cross section to NLO and above. There are, broadly speaking, two 
types of algorithm used for NLO calculations, differing in how they approximate the phase 
space and matrix elements in the neighborhood of divergent regions: 



1. The phase space slicing (PSS) method is based on approximating the matrix elements 
and the phase space integration measure in boundary regions of phase space so inte- 



gration may be carried out analytically 
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2. The subtraction method is based on adding and subtracting counterterms designed to 
approximate the real emission amplitudes in the phase space boundary regions on the 
one hand, and to be integra ble with respect to the momentum of an unresolved parton 



on the other 
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The phase space slicing method makes use of a combination of analytic and Monte Carlo 
integration methods, which has many advantages over a purely analytic calculation. The 
Monte Carlo approach allows one to calculate any number of observables simultaneously by 
simply histogramming the appropriate quantities. Furthermore, it is easy to tailor the Monte 
Carlo calculation to different experimental conditions, for example, detector acceptances, 
experimental cuts, and jet definitions. Also, with the Monte Carlo approach one can study 
the dependence of the cross section on the choice of scale and the size of higher order 
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quantum corrections in different regions of phase space. The basic challenge is to design a 
program which retains the versatility inherent in a Monte Carlo approach while ensuring 
that all of the required cancellations of singularities still take place. 

In this study, we use the phase space slicing method with one cutoff scale for which the 



universal crossing functions have been derived in Refs. |4l|, |42|, |43[ . The advantage of this 
method is that, after calculating the effective matrix elements with all the partons in the final 
state, we can use the generalized crossing property of the NLO matrix elements to calculate 
the corresponding matrix elements numerically without requiring any further effort. The 
validity of this method is due to the property that both the phase space and matrix element 
for the initial and final state collinear radiation processes can be simultaneously factorized. 
Below, we briefly review the general formalism for the NLO calculation in PSS method with 
one cutoff scale. 

The phase space slicing method with one cutoff scale introduces an unphysical parameter 
s m i n to separate the real emission correction phase space into two regions: 

1. the resolved region in which the amplitude has no divergences and can be integrated 
numerically by Monte Carlo method; 

2. unresolved region in which the amplitude contains all the soft and collinear divergences 
and can be integrated out analytically. 

It should be emphasized that the notion of resolved/unresolved partons is unrelated to the 
physical jet resolution criterium or to any other relevant physical scale. In the massless case, 
a convenient definition of the resolved region is given by the requirement 



Sij 



(Pi + Pj) 2 > Smin, (1) 



for all invariants = (pi + pj) 2 , where pi and pj are the 4- momenta of partons i and j, 
respectively. For the massive quarks, we follow the definition in Ref. [5l| to account for 
masses, but still use the terminology "resolved" and "unresolved" partons. In the regions 
with unresolved partons, soft and collinear approximations of the matrix elements, which 
hold exactly in the limit s m i n — *> 0, are used. The necessary integrations over the soft 
and collinear regions of phase space can then be carried out analytically in d = 4 — 2e 
space-time dimensions. One can thus isolate all the poles in e and perform the cancellation 
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of the IR singularities between the real and virtual contributions and absorb the leftover 
singularities into the parton structure functions in the factorization procedure. After the 
above procedure, one takes the limit e — > 0. The contribution from the sum of virtual 
and unresolved region corrections is finite but s m i n dependent. Since the parameter s min is 
introduced in the theoretical calculation for technical reasons only and is unrelated to any 
physical quantity, the sum of all contributions (virtual, unresolved and resolved corrections) 
must not depend on s m i n . The phase space slicing method is only valid in the limit that 
s m ; n is small enough that a given jet finding algorithm (or any infrared-safe observable) can 
be consistently defined even after including the experimental cuts. 

In general, the conventional calculation of the NLO differential cross section for a process 
with initial state hadrons Hi and H 2 can be written as 

d <% = E / dx 1 dx 2 f^(x 1 , t x F )f^(x 2 , t i F )da^°(x 1 x 2 ,ii R ), (2) 

a,b 

where a, b denote parton flavors and x±,x 2 are parton momentum fractions. f^(x,^p) 
is the usual NLO parton distribution function with the mass factorization scale y,p and 
da^ b LO (xi, x 2 , pLji) is the NLO hard scattering differential cross section with the renormal- 
ization scale The pictorial demonstration of Eq. (j2J) is shown in the upper part of 

Fig.m 

Contrary to the conventional calculation method, the PSS method with one cutoff scale 
will firstly cross the initial state partons into the final state, including the virtual corrections 
and unresolved real emission corrections. For example, to calculate the NLO QCD correc- 
tion to the H^-boson production via the Drell-Yan process, we first calculate the radiative 
corrections to W — > qq'(g), as shown in the lower part of the Fig. (TJ in which we split the 
phase space of the real emission corrections into the unresolved and resolved region. After we 
integrate out the unresolved phase space region, the net contribution of the virtual correc- 
tions and the real emission corrections in the unresolved phase space is finite but theoretical 
cutoff Smin dependent, which can be written as a form factor (denoted by the box in Fig. [I]) 
of the Born level vertex. 

Secondly, we take the already calculated effective matrix elements with all the partons 
in the final state and use the universal "crossing function", which is the generalization of 
the crossing property of the LO matrix elements to NLO, to calculate the corresponding 
matrix elements numerically. Once we cross the needed partons to the initial state, the 
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FIG. 1: Illustration of the PSS method with one cutoff scale to describe the processes with initial 
state massless quarks. Here, only half of the real emission diagrams is shown. In this paper, we 
will assign the particle's momentum such that the initial state particle's momentum is incoming 
to the vertex while the final state particle's momentum is outgoing. 

contributions from the unresolved collinear phase space regions are different from those 
with all the partons in the final state. These differences are included into the definition of 
the crossing function as well as the mass factorization effects, as shown in the middle part 
of Fig. [TJ Here, we only present the explicit expressions of the crossing function, while the 
definition and detailed derivation of the crossing function can be found in Ref. 42j. After 
applying the mass factorization in a particular scheme, the crossing functions for an initial 
state parton a, which participates in the hard scattering processes, can be written in the 
form: 



^scheme/ -i 
^ a \ i t^Fi °min) 



Nc 
2tt 



A a (x,[i F )\og 



+ £T heme (x,^) 



(3) 



where 



A a (x,fj, F ) = y] Ap^q (x, y, F ) 



(4) 
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5: chcme (*,M = 52Bgr(x,iM F ), (5) 
p 

and Nc denotes the number of colors. The sum runs over p = q,q, g. The functions A 
and B can be expressed as convolution integrals over the parton distribution functions and 



the explicit forms can be found in Ref. 42]. Although A a is scheme independent, B a does 
depend on the mass factorization scheme, and therefore so does the crossing function. 

After introducing the crossing function, we can write the NLO differential cross section 
in the PSS method with one cutoff scale as 

= / dx ^ dx ^fa 1 ( X l^F)f^ 2 (x 2 ,fiF)da^ b LO (x 1 ,X2,fiR) 

a,b J 

+ a 8 {fJi R ) [C^ 1 (x 1 ,n F )f^(x 2 ,n F ) + fa 1 (xi,^ F )C^ 2 (x 2 ,n F )] da^ b °(x 1 ,x 2 ). (6) 

Here da^ b LO consists of the finite effective all-partons-in-the-final-state matrix elements, in 
which partons a and b have simply been crossed to the initial state, i.e. in which their 
momenta — p a and -p b have been replaced by p a and p b , as shown in the Fig. [TJ The 
difference between da^ b LO and da^ b LO has been absorbed into the finite, universal crossing 
function C^(x,fi F ). Defining a "effective" NLO parton distribution function J-^{x) as 

= (x,im f ) +a s (fx R )C» (x,im f ) +0(a 2 s ), (7) 

we can rewrite Eq. ([6]) in a simple form as 

= E / dx l dx 2 ^{x l )T b H \x 2 )dal LO {x l: x 2 ). (8) 

a,b •* 

III. TREE LEVEL MATRIX ELEMENTS 

The matrix elements of the scattering process gb — > tW can be written in terms of the 
following twenty standard matrix elements that contain the information about the Dirac 
matrix structure: 

M LO = M$ ) + M%, (9) 
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FIG. 2: Tree-level Feynman diagram of the gb — ► iPF process. 



where 



(0 

LO 



+fePt,^u + fjPb^Pb,v + flPb,nPt,v + flPt,flPb,v + fl Pt^Pt,u 



+fi&Pt,iM7u + fl 7 Pb,»Pb,v + fl 8 Pb^Pt,u + f\ 9 Pt,^Pb,u + fmPt,tfH,v) 



where i = s(t). Here, g s (g) denotes the coupling strength of the strong (weak) interaction, 
respectively, T" 1 - are the color matrices in the fundamental representation, and pi = (1 — 
7s)/2 denotes the left-handed projector. The tree-level Feynman diagrams of the gb — > tW 
process are shown in Fig. [2j The non-zero form factors of the tree-level matrix elements are 
given by 

2 

s 



rs rs r s rs r> r s 

J3 ~ J A ~ J5 ~ J 11 ~ z il2 

2 



(11) 



fl = V 



The kinematics variables used in this paper are defined as follows: 

s = (Pg+Pb) 2 = 2pg-Pb, 
t = (Pg-Ptf = h + rrfi, 
u = (Pg- Pw) 2 = ui + mw, 



(12) 

(13) 
(14) 
(15) 



where t\ = —2p g ■ p t , u\ = — 2p g ■ p w , and s + ti + Ui = 0. 



IV. VIRTUAL CORRECTIONS AND RENORMALIZATION 



Now let us calculate the virtual corrections to the tW associated production. Fig. [3] 
shows the Feynman diagrams of the one-loop QCD virtual corrections, where Vi denote the 
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FIG. 3: Feynman diagrams of the one-loop virtual corrections where the red lines denote the top 
quark. 

triangle loop corrections, B{ the box loop corrections, while Si the bubble loop corrections. 
In this work we adapt the on-shell renormalization scheme such that the external self- 
energy corrections are cancelled by a set of counterterms which will be discussed later in 
this sections. 

At the NLO the relevant one-loop virtual corrections contribute only through their in- 
terference with the lowest-order Born amplitudes. The interference, i.e. 9ft (-Mo-M. virtual J > 
gives rise to the order of a 2 s contributions. One can further categorize the virtual corrections 
into three classes, according to their color structures depicted in Fig. HJ (a) V±, V{, B 1: B 2 ; 
(b) V 3 , V^, Si, S2; (c) V 2 , Vi, B 3 . The color factors (CF) of those three classes, normalized 
to the Born level matrix element square, are given by 



CF a 



1 



CF h 



Nc 
2 



Nc 
2 



(16) 



2iV 2 2N C 

The same color structures also apply to the real radiation corrections. Below we will make 
use of the color order to organize our calculation. 

To calculate the virtual corrections, we follow the Passarino-Veltman procedure 
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53|: 
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(a) (6) (c) 

FIG. 4: Color factor of the interference of tree- level and one-loop matrix element. 

one can reduce any set of vectorial of tensorial one-loop integrals to a few scalar functions. 
However, since the loop corrections exhibit not only the UV divergence but also the IR 
divergence, either soft or collinear or both, one needs to work out those scalar function 
analytically. Recently, a basis set of infra-red and/or collinear ly divergent scalar oneloop 
integrals is constructed, and analytic formulas, for tadpole, bubble, triangle and box inte- 



grals, are given explicitly in Ref. [54|. The scalar functions used in this work can be found 
in that paper and those results are confirmed by our hand calculation. Our final results 
of the virtual corrections are very complicated, but the divergent pieces can be factorized 
out in a simple form. We present the full expressions and the divergent poles of the form 
factors in the Appendix O and [Dj respectively. We further distinguish between the UV and 
IR divergences in order to keep track on the renormalization and factorization. 

Consider now the UV divergence first. We define Auv(-M-i) to be the UV pole part of 
the corresponding amplitude M. i . Using this notation, we find 

A UV (M Vl + M V () = %C. (--L) ±- Muh (17) 

A vv (My, + M vi ) = ^C, (^f) —M LO , (18) 

A UV (M V , + M VS ) = £ C . (f - ±. Mlo , (19) 

Auv (M St + Ms,) = (f - ji-) Am 10 + A Sl , (20) 

where C e = (4nfi 2 /mf) 2 T (1 + e). Note that As 1 in Eq. [201 represents the non-factorization 
loop corrections originated from the Feynman diagram Si, see Eqs. ( 1C861IC88I) . The UV 
divergent poles of ,4^ are given in terms of the form factors as follows: 



ff = =« J-, = 2/l f ■ = -L - i ) .21. 

tx euv e uv V h l \ 
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As to be shown later, these non-factorisable divergences will be exactly cancelled by the top 
quark mass renormalization. 

The renormalization is preformed in the MS scheme with the top quark mass defined on 
shell. As required by renormalization group arguments, the renormalization of the funda- 
mental propagators and interaction vertices of the theory reduces to introducing countert- 
erms for the external field wave functions of top quarks and gluons (SZ t , SZ G ), for the top 
mass (5m t ), and for the strong coupling constant (5g s ). We renormalize the fields of the 
gluons, G a ,/i) of the bottom quark, ip b , and of the top quark, ip t , all in the on-shell scheme, 
i.e. the wave-function renormalization constants SZc,b,t, defined by the transformations 



(22) 



are adjusted to cancel the external self-energy corrections exactly. Distinguishing between 
divergences of UV and IR origin, these constants can be written as 

a, „ fN c 1 \ / 1 2 



SZ t = -—C e — - — ) — + — + 4 , (24) 
4tt V 2 2N C ) \e uv e IR 



a. 



where rif = 5 is the number of light quark flavours 

n i n f 1 in I \vivn f / » i ^ / n i -n ■ L- iiimcl' m i w 1 fl in I ivirn ct-v/uwi /■/ mnlnwi 1 ■.ic 1 > ; ■ t, \ ■■ \ \ -. \ 



Denoting the bare top quark mass and the bare strong coupling as m° and a®, respectively, 



we introduce the renormalization parameters in the transformations 

m° t = m t + 5m t , a® = a s + 5a s . (26) 

We defined the subtraction condition for the top quark mass m t in such a way that m t is 
the pole mass, in which case the top mass counterterm is given by 

5m t a s ( Nc 1 \ / 3 



mt ATi Ct \2 2N c )\e uv +4 )' (2?) 

Finally, for the renormalization of a s , we use the MS scheme, modified to decouple the top 
quark from the running of the strong coupling a s (//) . It gives rise to 

Sa -- a ' Wr(i + e )(^-i^-L + Hc. 1 



47T ^ \ 3 3 / €uv 3 euv ) 



11 



® 8m t 
-<WW\, 



6a, 



*VWW 



*wwv 



sz b 



*vww 



FIG. 5: Depicted Feynman diagrams of the counterterm contributions. 

where the first term inside the braces originates from light quark and gluon loops, which 
is independent of fi, while the second term (proportional to C e ) originates from the top 
quark loop in the gluon self-energy that is subtracted at zero-momentum transfer. In such 
a way the running of the coupling a s (//) is generated solely by the finite contributions of 
the light quark and gluon loops, while the top quark contribution is absorbed completely in 
the renormalization condition and thus decouple effectively. 

The contributions of the counterterms and the external self-energy to the one-loop matrix 
element are depicted in Fig. [5j The matrix elements is given by 



(^T + \ 5Zb + \ 5Zt + \ 5Zg I Ml ° + M ^ mt ^ 



(28) 



The first term accounts for the external self-energy corrections as well as the counterterms 
of the strong coupling renormalization, while the second term proportional to 6m t comes 
from the renormalization of the top quark mass. After insertion of the top quark mass 
counterterm, see Fig. [5j we obtain the form factor of M. (6m t ) as follows (where the color 
factor y — is not shown explicitly): 



fSmt 
J2 



fSmt " "t 

^6 — 7T 



9 

mi 



' 12 



+ 4 



+ 16 



12 



Sm t 



(29) 
(30) 
(31) 



In order to check the renormalization, we will combine M. (Smt) with those bubble, 
triangle and box loop corrections, therefore the renormalized matrix elements are given 
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by 



Mren = M virt + M C T, (32) 
"23 

M , = Ms i + ( Mv < + M V + M Bi ) + M (5m t ) , (33) 



=1 i=\ 



Mct = ( ^ + \8Z b + l -5Z t + \bZ G ) M LO - (34) 



It is straightforward to check the UV divergences indeed cancel out, but the IR (soft and 
collinear) divergences still remain to be cancelled by the real radiation corrections. The IR 
poles of the virtual corrections are given by 

a wirt a a a s ( 13 3 . -h 3 s 1 m\-m 2 w -u 
A v /r = Mlo x — — =- + ln^ + hi—o- o ln — T 

10 1 /11JV C 2 
3e/H 2e/R \ 3 



n/H. (35) 



V. INFRA-RED SINGULARITIES IN THE REAL EMISSION CORRECTIONS 

In the one cutoff phase space slicing method, the theoretical cutoff parameter (s m j n ) 
is introduced in order to isolate soft and collinear singularities associated with real gluon 
emission sub-processes by partitioning the phase space into soft, collinear and hard regions 
such that 

|A*f = l-^toft + l-Mtollincar + l^lLd ■ (36) 

In the soft and collinear regions the cross section is proportional to the Born-level cross sec- 
tion. Using dimensional regularization, we can evaluate the real gluon emission diagrams in 
n-dimensions under the soft gluon approximation in the soft region, or the collinear approx- 
imation in the collinear region, and can integrate out the corresponding phase space volume 
analytically. The resulting divergences are cancelled by virtual corrections or absorbed into 
the perturbative parton distribution functions in the factorization procedure. 

In this method, a pair of partons with momenta pi and pj is defined to be unresolved if 

\2pi -Pj\ < Smin, (37) 

with s m in small compared to the hard scale of the process. This condition can occur if 
either pi and pj are collinear, or if one of the two is soft. When the scattering amplitude 
involves a complicated color structure, one needs to decompose the scattering amplitude 
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FIG. 6: Feynman diagrams of process W — > i&g at Born level. 

into color-ordered sub-amplitudes. Below we will use the tW associated production process 
to illustrate this point. 

To calculate the real emission corrections, we first cross all the initial state partons to 
the final state and the final state W boson into the initial state, i.e. at the tree level, 

gb w - t w + ti g . 

The tree-level Feynman diagrams after crossing are shown in Fig. [6j which give rise to the 
following matrix element 
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M = -t-^=g s T a e»(Wy*(g)u(t) 



p z — m t q z — mi 



v(b), (38) 



where q = p g + pi and p = Pt + P g - The cross symbol "x" in the figure indicates the 
possible places where the additional parton can be radiated from. There are two types of 
real radiation corrections: 

(1) W + tbgg; (2) W + tbqq. 



The former exhibits both soft and collinear divergences, but the latter can only have the 
collinear divergence. 



A. Color ordered amplitude 

For a systematic extraction of the infrared singularities within the one cutoff method, we 
organize the amplitude in terms of color-ordered sub- amplitudes. In Fig. [71 we present all 
the real emission diagrams which give rise to the final state of tbgg. The color coefficients 
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FIG. 7: Feynman diagrams contributing to tbgg final state. 



of the diagrams are given as follows: 



a,b,c oc T a T\ 



iarpb 
-ibrria 



d,e,f oc T b T a , 
g,h oc T a T b and T 6 T a , 



where we have applied the identity 

^ jabcrpc rparpb rpbrpC 

in the diagrams (g) and (h). Thus, we can decompose the amplitude Ai(W + 
following, 

M(W + -> tbgg) = T a T h M \ + T b T a M 2 , 

where 



.Mi = .M {a) + M {b) + .M (c) - M {9) - M {h \ 
M 2 = M {d) + M {e) + .M (/) + M {9) + .M w . 
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The color indices are assigned as shown in Fig. The momentum q i} qij, p, and Pij are 
defined as 

Pi = Pt + P 9i , Pji =Pt+ P 9l + P 9] , 
qi=Pb + P 9i , qji =Pb + P 9i + Pg 3 ■ 



The amplitudes of diagrams (a, b, c) are given as follows (with a color coefficient T a T b )- 



= l -^= g y(w + y*( gi y*(g 2 ) 



x u(t) 



p - 4\2 + m h - q , 2 + m b 

111 L o 1 2 2 ^ a 

9ia - m b %- m b 



= z^L g y(w + y*( gi y*(g 2 ) 



x u(t) 



j>\ + rn t - j 2 + m b 

Tp 2 2 T/i-* £ 2 2 ^ 

p{ -mi * qi-ml 



m® = t^= 9 y(w + y*( gi y*(g 2 ) 



x u(t) 



^i + ra t ^2i + rat p 

^P-2 ™2^°"„2 ™2^M ^ 



0. 



(6). 



pj - m£ P21 - ra t 

The amplitudes of diagrams (d, e, f) are given as follows (with a color coefficient T b T a ^ 



(41) 



(42) 



(43) 



= z^= g 2 y(w + y*( gi y*(g 2 ) 



x -u(t) 



- (/21 + ra b - j/i + m b 

'Yu-LL 2 2 2 2 7p 

521 - m fe 51 - ra^ 



x -u(t) 



, v ^2 + rn t - 4i + m b 

^ u -2 ™2^P £ „2 ra 2 ^ 



P2 - raj 



5f - ra^ 



(5) 



x u(t) 



^2 + ra t ^12 + ra t 

^ /fT --2 ^,2^P^2 ^2^P- L 



(44) 



(45) 



(46) 



P2 - ra? pij, - mt 

Finally, the amplitudes of diagrams (g, h) are given as follows (with a color coefficient 
-if abc T c ): 



= i^=g 2 y{w + y\ gi y*{g 2 ) 



x u{t) 



- 4k + m b 

lp,Ph— 2 2~ 7a 

Tk ~ m b 



1 



(47) 
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= i^= g y(w + y*( gi y*(g 2 ) 



x u(t) 



Pk +m t 

la— 2^^ F L 

. P k - m t 



(b) -V™, 



where 



The squared amplitude, after summing over colors and spins, are given by 



J2\M(W + ^fbgg 



color 
spin 

N 2 C -IN C 



(48) 



(49) 



(50) 



where summing over the final state degrees of freedom is understood in |.Mj| 2 . The factor- 
ization of soft and collinear singularities for color ordered amplitudes has been discussed in 
the literature mainly for the leading color terms (O(N)) 4l|, |42|. For our calculation of the 
one cutoff phase space slicing method, we will have to extend these results to the sub-leading 
color terms (0(1/ N)). 



B. Soft singularities 

The soft gluon behavior for an color ordered sub-amplitude is very similar to the soft 
photon behavior of QED amplitudes. In QED, the soft photon couples to a charged fermion 
line, resulting in an eikonal factor multiplying the hard process. The key point is that only 
the Feynman diagrams with soft photon coupled to the external charged fermion lines will 
contribute in the soft photon limit. Technically, no matter how many photons are radiated 
out from the charged fermion dipole, the soft eikonal factor only knows about the external 
momentum since there is no photon-photon interaction in QED theory. In QCD, the soft 
pattern becomes much complicated due to the non-Abelian interaction. However, it has 
been shown that the color ordered sub-amplitudes do exhibit a factorization of the soft 



gluon singularities as in the QED [55l . [56|. This is because the partons are ordered and 
form well defined color charge lines to which the soft gluon can couple. The soft gluon 
behavior depends only on the momenta of the external color charged lines to which the soft 
gluon couples, and is independent of the number and type of other partons in the process. 
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FIG. 8: Momentum label of the color ordered amplitudes. 



Similarly, the soft factor is independent of whether or not any color singlet particles such as 
electroweak bosons are participating in the hard process. 

1. Factorization of the color ordered amplitudes in the soft limit 

In the limit of one gluon being soft, each term in the squared amplitudes (cf. Eq. l50l) can 
be written as a factor multiplying the Born level amplitude square. Since t and b are the 
external fermions which cannot be soft, there are only two sources of the soft singularities: 
either from gi or from g 2 . When both of them are soft, it goes beyond the NLO (0(a 2 )) 
and will not be considered in our calculation. 

First consider | i, 2 1 2 - The momentum configurations which respect the color order are 
shown in Fig. [HI In the soft gluon limit, |.Mi| 2 can be factorized as follows: 

N c~ 1 N c , u ,2 .i-o g 2 s N c 2 

— 2 2~\Mi\ * —^-f(t, 91,92) \Mo\ , (51) 

^ g -^f(9i,92,b)\M \\ (52) 

where the eikonal factor f(a, s, b) is defined as 

ft h\ 4Sab 4m * 4m b w 
f(a,s,b) = 2 — . (53) 

SasSsb S as S sb 

Here we define Sjj = 2pi ■ pj for both massive and massless partons. The color factor 
(Nq — l)/2 has been absorbed into the Born level amplitude square. Similarly, |^W 2 1 2 can 
be factorized in the soft gluon limit as follows: 

2 A T 

► — 2-f(t,g 2 , 91) \Mq\ • (55) 
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2 

Now let us consider the QED-like term \A4i + A4 2 \ ■ This interference term exhibits the 
canonical eikonal factorization in the soft limit. The soft gluon only knows about the external 
color charged fermion lines, leading to the following factorization 

-^2^' + ^| 2 ^ -4-J(^)\M^, (56) 

■ !! /(«,»,!>) |A1„| 2 . (57) 



2N C 



2. soft singularities of \M.\\ 2 



In order to slice the phase space we introduce the following Heaviside functions 0(i, j, k) 
and k): 

k) = Q(sij + s jk - 2s min ), (58) 
k) = Q{2s min - Sij - sjk). (59) 

With the help of 0-functions we can split the phase space of the real emission corrections 
into four regions: 

1 = (6*9192 + 6*9192) (@9192& + ©9192&) 

= ©*9l92©9i9 2 5 "I" ®^l92©9i9 2 S + ©*9i92©9i9 2 5 — @*9i92©9i9 2 5- (60) 

The first term denotes the region where two gluons are both hard (no soft singularities), 
but one should keep in mind that there might exist the collinear divergences, depending on 
the masses of external fermions. The second term denotes the region where the gluon g 2 
is soft but the gluon g x is hard. The third term denotes the region where g 1 is soft and 
g 2 is hard. The last term denotes the region where g 1 and g 2 are both soft. As discussed 
before, this region only contributes when one calculates the quantum corrections beyond the 
next-to-leading order and therefore is ignored in this work. The negative sign ahead of the 
fourth term is to get rid of the double counting from the regions described in the second 
and third terms. 

When gi is soft, the eikonal factor f(t,gi,g 2 ), 

f{t, gi ,g 2 ) = — — y - -3-, (61) 
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can be calculated in the center frame of (t, g?) system where we choose the explicit momen- 
tum of final state particles as follows: 

p t = (E t ,0,0,p t E t ), 
Pg2 = (p t E t , 0, 0, -p t E t ) , 
p gi = (E gi ,0,E gi sin9,E gi cos 9) . 

Here, E t , E gi and E g , 2 is the energy of the top quark, the soft gluon [g{] and the hard gluon 
(g 2 ), respectively, and 9 is the angle between the soft gluon g 1 and the top quark. It is easy 
to show that the eikonal factor can be written as 



1 



El E t 



where 



mi 



1 + cos 0)(1 -Acos 9) E t (l - p t cos 9) 2 



(62) 



s tg2 = 2/^ 2 (l + /3), 

Stg 2 



Pt 



s tg2 + 2m 2 ' 
s tg2 + 2m 2 



2\/ s tg2 + 



mi 



Now let us calculate the phase space boundary condition with the choice of momentum 
above. The sum of s tgi and s gig2 is given by 

s t91 + s gig2 = 2E t E gi (l + p t ) = 2E gi ^s tg2 + ml (63) 

Substituting Eq. [63] into the soft phase space boundary condition, we obtain the upper limit 



of E. 



91 • 



E 91 < 



V Stg 2 + 



m 



— rpmax 
2 ~ 91 ' 



In the soft limit, the phase space is also factorizes as 

1 (4vr 



'87r 2 r(l-e) 91 
Then the soft gluon contribution is given by 



Elr 2t sin 1 -^ 9 dE n , d9. 



(64) 



I sof %9i,92 

~9 2 S N< 



£9192 



c 



{^ d y f(t, 91,92) 



d 3 9l 



{27r) d - 1 2E, 



31 



N, * s 



le- 



gs 1 (47T/X 



2\ e 



2 87r 2 rd-e) 



run <r.r 
E S1 



^i 1_2e d^i 



x / sin 
'0 



l-2e 



,1 

Et 



mi 



mi 



1 + cos 9) (1 - p t cos 9) E t {l-p t cos9Y 



c\9. 



(65) 
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Using the formula listed in Appendix |A] we evaluate the integral in Eq. [65] and obtain the 
soft factor I so ^ (t,gi,g 2 ) as following, 



I^\t, gi ,g 2 ) = N ( 



l 



c 



167r 2 rd -e) \ s 



47T/X 



2\ e 

d 



s tg2 + m t 



ln( l + ] +21n2-l 

mi 



7T~ 

~6 



+ 21n 2 2-21n2 + 



21n2 + 



s t92 + 2m 2 



In 1 + 



2 

mf 



--hr ( I . + ] -2Li, 



2 \ mi J ' \s tg2 +mi 
Similarly, one can calculate the soft singularity when g 2 is soft, yielding 



(66) 



i s °f\ 9l) 92 ,b) = Je 



N, 



9i92b 

9l 



9 2 S N C f 2Y f( t, d 3 g 2 
~rr~ {Vd) f {9i, 92,0); 



(2vr)^2£ S2 



c 



16vr 2 r(l -e) Vs 



2 41n2 , 2 7i 
-a + 41n 2 2- - 



3 mm 



mm / \ ^g\b 
2 



~IR 



(67) 



This result is for = only. If m;, ^ 0, then the soft factor can be obtained from Eq. 



J*%i,<fc,&) =N< 



o 



i6tt 2 rri -e) V s 



2\ e 

d 



s gi b + mi 



1 1 



e 2 e 

7T 2 



In 1 + 



S 9lb 

2 

mi 



+ 2 In 2 - 1 



21n 2 2-21n2 + 



21n2 + 



s gi ~ b + 2m 2 b 



'gib 



In 1 + 



'gib 

2 

mt 



- In 2 1 1 + 8 » b 



mt 



-2Li 2 



'gib 



s gi b + mi 



(68) 



3. soft singularities of \A4$ 



The phase space can be splitted into four regions: 



1 — (@t329i + @*929i) (©92916 + ®<?2<?ib) 

= ©*929i©929i5 + ®tg 2 gi®g 2 g 1 b + ©*929i©929i6 — ©*929i ©9291&" 



(69) 



The soft singularities of |_Axf 2 1 2 can be derived from the results of |-Mi| 2 by making the 



substitution g x <->• g 2 . 
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4- soft singularities of \M.\ + 



The phase space can be splitted into 

1 = (e tfll6 + e tfll5 ) (e tg2 - b + e tg2l ) 

= ®tgib®tg 2 b + ®tgib®tg 2 b + ®tg 1 b®tg 2 b ~ ®tg 1 b®tg 2 b- 



(70) 



The first term corresponds to the two hard gluons, the second term denotes E gi — > 0, the 
third term denotes E g2 — > and the fourth term denotes both E gi — > and E g2 — > 0. Since 
the eikonal factor is only determined by the two external fermion lines (t and b), the soft 
factors of g 1 and g 2 should be the same, i.e. 

I s °f t (t,g 1 ,b) = I soft (t,g 2 ,b). (71) 

The soft factor can easily be derived from that of &t gi g 2 l-^i| 2 by making the substitution 
(?2 — > b, because the soft gluon contributions only know about the kinematics. It gives rise 
to the soft factor as following, 



r oft (t, 9l ,b) = i soft (t,gi,92) 



92^b 



c 



9t 



1 



X 



i67r 2 rd 
1 



47T/i 



2\ e 
d 



Sfb + m 



Z IR 



s tb 



In [ 1 + ] +21n2-l 



■— + 21n 2 2-21n2 + 
6 



21n2 + 



Sth + 2m 



Stb 



In [ 1 + 

mi 



•-In 2 



1 + ?* ] _ 2Li 2 

ml 



Hb 



so, + mi 



(72) 



C. Collinear singularities 

As mentioned above, there still exists collinear divergence even for a hard gluon radiation. 
Below we will further slice the collinear region of the phase space to derive the collinear 
divergences. 



1. Collinear singularities of \M.\\ 

We further slice the hard phase space of \M\\ as 

©tfllfl2© S i S2 6 l-Ml | 2 = { < dtg l92 Q gig2 b - C Ml + C Ml ) \M\\ 2 , 
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(73) 



where 



~l~ 2s m j n ) ( y Sg 1 g 2 2s m i n ) {Smin S 3 2 b) " (74) 

The first term in (7^ represents the phase space region that neither g 1 or g 2 is not soft but 
gi || <7 2) while the second term denotes the phase space region of g 2 \\ b. 

In the collinear region of g± \\ g 2 , the color ordered matrix element square exhibits the 
following factorization property, 

^^_± {Mif »^ 9jNc f99 ^ 9 {Mof (?5) 

where 

^(q-J- ' + f + ('-«' -*!*, (76) 



S 



Si 92 



with the kinematics defined as 

Ph = p gi + P g2 , p gi = &h, P 92 = (1 - Pfc- (77) 

The phase space integration in the collinear region can also be written in the following 
factorization form 

d n p 9l d n p g2 



(2 7 r) n - i 2E 9l (27T)"- 1 2E, 



92 



1 (M) [wa-or^x-^-. (78) 



"" 167r2r(l- e ) ™ 2SV S7J ^ 1S2 s (2vr) d - 1 2^' 

Bearing in mind that the integration is limited by the conditions s gig2 < s min , s tgi > 2s min 
and s g2 i > 2s min , we integrate the matrix element over the collinear phase space and obtain 
the following collinear factor, 



= a s N c (4xnY[ 2 
Ml ~ 47rr(l-e) \s mm ) \e IR 



2 ' ^ 2s ■ 2s • II 
In h In 1 



2??7.£ 27r 67 2 2s m j ra 2 



_ ^2 _ ^2 ^^mn , (?g) 

s t<n 3 9 s tfll s gib ) 

In the collinear region of g 2 \\ b, the color ordered matrix element square exhibits the 
following factorization property, 

l-Mil 2 ^ 9 ^r-« |^o| 2 , (80) 
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where 



2 l + e-e(l-ZY 



(81) 



Sg 2 b (1 - 

Again, after integrating over the collinear phase space, we obtain the following collinear 
factor, 



4tiT(1 -e) Vs. 
2. Collinear singularities of \J\A2\ 2 



(D 2 < 


\1_ 


2s m j n 3 
In h - 






- s <n5 4 J 



7T 2 7 , . 

1 In' 



2 s 



2 ^°min 



'gib 



(82) 



1 1 2 1 1 2 

Similar to the case of \Mi\ , we split the hard phase space of \M.2\ as 



®tg 2 gi®g 2 g 1 b 



®tg 2 g 1 @g 2 g 1 b ~ Cm 2 + Ca4 2 



l-M 2 l 2 , 



where 



C*x 2 — (st 92 — 2s m i n ) — 2s m inj (s 



~l~ ® (^t& 2s m j n ) (.Sg 2 gi 2s m j n ) (s m j n ^gib) ■ 



(83) 



(84) 



The first (second) term in Cm 2 represents the collinear region of #2 || 9\ (gi \\ b), respectively. 

2 

The collinear factors can be derived on the analogy of those of \Mi\ , 

U M 2 - U Mi ' ^M 2 - U M 2 ■ y 6b ) 
3. Collinear singularities of \M.\ + M.2\ 2 



We split the hard phase space of \Mi + M 2 \ = M-i+2 as 



@i<7i6@i(;2b l^l+2| 



®t gi b®tg 2 b ~ Cmi +2 + Ca1i+2 



I At 



1+2 1 



where 



C-Mi+2 — ® 2s m j n ) (s ff2 ^ 2s m j n ) (s m j n 

"I - {ptg 2 2s m j n ) (s^fe 2s m j n ) (s m j n ^^s) • 



(86) 



(87) 



The first (second) term in Cm 1+2 represents the collinear region of g\ || b (g 2 \\ b), respectively. 

In both collinear regions, the matrix element square can be factorized as 

,2 

2 3i||0 i 

1 1+2 1 



N C ~ 1 1 I « A 



® s tigi^g 1 a a 1 2 



2iV c 



2iV c 



(88) 
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Integration over the collinear region gives rise to the following collinear factor, 



U Ml+2 " 

a* 1 



C<qgi-^q 

1 



4nN c T(l -e) Is 



2\ 2 



€lR 



In h - 



2mjf 7r 2 7 2 2s m j n 



(89) 



^. Collinear singularities of W — > t&gg 



Now let us consider the last piece of the real emission correction, VT + — > t&gg, which only 
involves the collinear divergence. We split the phase space into two parts, 



(90) 



where the first (second) term represents the finite (collinear) region, respectively. In the 
collinear region q \\ q, the matrix element square can factorized as 



where 



\M (tbqq)\ 



f qq->9 (0 



2 q\\q g ( 



f q ^ 9 \M \ 2 , 



2 e + a-o 2 -£ 



1-6 

Integration over the collinear region gives rise to the following collinear factor 

1 /47m 2 Y f 2 10' 



(jqq-^g 



47rr(l-e) \s r , 
where rif denotes the light quark flavors. 



x 



n f , 



(91) 



(92) 



(93) 



D. complete result of the IR singularities of the real emission corrections 

The complete soft and collinear singularities of the real emission corrections to the process 
gb — > tW are given by 



vir = (Is + Ic) l-Mol 2 , 



(94) 



where the soft factor I s is given by 



J * = 2 



I s °f l (t, 9l , b) + 1^ (t, g 2 , b) + (t, g 2 , 9l ) 
{g 2 , 9l , b) + J"* (t, g u b) + 1^ (t, g 2 , b) 



(95) 
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while the collinear factor Ic is given by 
Ic = 



r^gg^g I r^qg^q 



r<99^q i r<qg-*q , r"i9i^q , r<i92^q 



(96) 



The factor "2" in above equations accounts for the two identical gluons in the final state. 
After summing over both soft and collinear factors, we obtain the IR divergences of the real 
emission corrections as following 



&IR 



a* 
Air' 



26 
Wr 



6 



In 



em V 3 



Stg_ 

2 

mi 



mi 



£ir 



mi 
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+ finite > | .M 1 



(97) 



It is clear that the IR divergences in the real emission corrections exactly cancel the UV 
divergences in the virtual corrections after we cross the gluon and bottom quark into the 
initial state and the W-boson to the final state, i.e. 



-u 



s, 



m\ — m^ — u. 



(98) 



VI. RESOLVED REAL EMISSION CORRECTIONS 

The real emission corrections in the resolved phase phase region are finite and thus can be 
calculated numerically in 4 dimensions, using the canonical Monte Carlo method. With the 
implement of the phase space slicing conditions described above, the integration will depend 
upon the input value of s m i n . Since the cutoff s m i n is introduced in the calculation only 
for a technical reason and is unrelated to any physical quantity, the inclusive rate must not 
depend on it. In other words, the sum of all contributions, virtual, resolved, and unresolved 
corrections must be independent of s m i n . This is the case as long as s m i n is small enough so 
that the soft and collinear approximations are valid. However, numerical cancellation in the 
Monte Carlo integration becomes unstable if s m i n is too small. Furthermore, the jet-finding 
algorithm and other infrared-safe experimental observables should also be defined in a way 
such that they are consistent with the choice of s m i n . In practice, one wants to choose the 
largest s m i n possible within these constraints in order to minimize the processing time of the 
Monte Carlo integration program. Since the phase space slicing conditions highly depend on 
the kinematics of the final state particles, it is better to calculate the resolved real emission 
corrections term by term in Eq. [50] to achieve better cancellation of the s min dependence. 

26 



VII. CONCLUSION 



In this paper we present a detailed calculation of the NLO QCD corrections to the tW 
associated production using the one cutoff phase space slicing method. The corrections to 
the this process have been calculated independently by two groups 
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13511 , a number of dif- 



ferences were found. Therefore, a third party calculation is needed to make a crossing check. 
Furthermore, the fully analytic expression is still missing in the literature. In this note, we 
first calculated the virtual corrections using the dimensional regularization. The final re- 
sults are shown completely in terms of the Passarino-Veltman scalar functions. We then use 
the one cutoff phase space slicing method to calculate the infrared singularities of the real 
emission corrections. The soft and collinear factorizations of the color ordered amplitudes 
are shown in details. The phenomenology study will be presented in the forthcoming paper. 
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APPENDIX A: USEFUL FORMULA 

The integration over the soft and collinear regions can be easily evaluated. Below we list 
the useful formula in our calculation [57]. All the integrals we need are of the form 

rw> = f\w sinl ' 2£gl f ha sin " 2e6>2 fAn 

J 1 (o + 6cos0i) fc 7 o 2 (A + Bcose 1 + Csm9 1 cos9 2 ) 1 ' 1 ; 

The first point to notice is that there are four classes of integrals, depending on the collinear 
structure. The [ab] and [AB] variable can both be either "collinear" or not, which yields four 
combinations. In the case of the [ab] variable, collinear divergence appears when a 2 = b 2 so 
that (a + 6cos#i) — > a(l ± cos#i). Then the zero occurs at the edge of the 8\ integration 
region and is not integrable. For the [yl.BC] variable the same comments apply for 82 when 
A 2 = B 2 + C 2 . For our needs, the soft eikonal factor does not involve the angle 82, therefore 
the integration becomes more simpler. 
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When A 2 ^ B 2 + C 2 , and b = —a, we use 



7T 



a(A + S) 



h In 

e 



(A + B) 1 



A 2 - B 2 - C 2 



A-VB^+C 2 
A + B 



In 2 



A + VB 2 + C 2 



+2Li2 |-^Sl-2Li 2 



A- VB 2 + C 2 ^ 

B - VB 2 + C 2 
A + B 



While when b ^ —a we use 



j(o,i) 



7T 



Vfi 2 + c 2 



In 



A + VW+C 2 



+2e 



/(0,2) 



Li 2 



2vr 



A-y/WTC 2 
1 



2V5 2 + C 2 



A + VW+C 2 



A + VW+C 2 



A-VW+C 2 



A 2 -B 2 - C 2 



1 + e 



A , / A + V5 2 + C 2 
: In 



V£ 2 + c 2 XA-tJW+C 2 



Finally, when A = B 2 + C 2 , and 6 = —a, we have 



/(i.D = _i_L fA + B 



eaA\ 2A 



-1-e 



1 + e 2 Li 2 



A-B 
2A 



+ 0(e 2 ), (A2) 



+ 0(e 2 ), (A3) 



(A4) 



(A5) 



APPENDIX B: PASSIANO-VELTMAN SCALAR FUNCTIONS 



Here we list out the scalar functions used in this calculation. The general scalar functions 
are defined as follows: 



A (m 2 ) = fj, 



2c 



d n q 1 
(27r) n q 2 — m 2 
d n q 



B (pl;m 2 v m 2 2 ) = fi 2e J 
Co {pl,pl,pl;m 2 l} ml,ml) 



(27r)™ [q2 - m 2 ] [(g + pi ) 2 - m 2 ] 



1' 



2c 



d n q 



(27T) n [g2 _ m 2j + p ^ _ m 2j ^ g +pi+ p2 f _ m 2j 



n / 2 2 2 2 2 
A) (Pl,P 2) P3^4; s 12,S23;^l 



2 2 2 2 

,m 2 ,m 3 ,m 4 



) 



(Bl) 
(B2) 

(B3) 
(B4) 



2c 



d n q 



(2vr) n [g2 _ m 2j j-^ _|_ ^ _ m |j [(g + pi + t? 2 ) - m§] [(<? + pi + p 2 + P3) - m 2 ] ' 

where n = 4 — 2e, /i is the scale introduced so that the integrals preserve their natural 
dimensions. We follow the notation of Ref. j^] to calculate the scalar integrals in the 
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spacelike region, p 2 < and s^- = (p* + £>j) 2 < 0. The analytic continuation is performed by 
restoring the ie, 



Pi 



Pi + ie, Si 



Sij + ie, m i 



to,- — ie. 



The divergent three-point scalar functions used in this work are given as follows: 
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where 
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r(i + e). 



The divergent four-point scalar functions are given as follows: 
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it 



m: 



-2Li 2 



mr — m 



w 



m\ — t 



2Li 2 



mr — m 



w 



mf — u 



If 



(0,0,m 2 t ,m 2 w ;s,u;0,m 2 t ,0,0 



1 



1 



167T 2 £ s(«-m^) [2e 2 e 



3 1 



2 In 1 - 



u 



mt 



In 



m: 



-In 1- 



rn 



w 



mi 



-2Li 



u - r?% 

2 | 2" 

u — mf 



TO 2 



f21n(4)lnfl-^ 



TO 2 



In 2 I 1 



TO 2 



3* 



D B 2 



(B8) 



(B9) 
(B10) 



The scalar functions exhibit certain symmetries under interchange their arguments (either 
rotation or inversion). The symmetry properties of the triangle- and box-loop integrals used 
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in this calculation are listed as follows 54| : 
(1) Cq- function 



C (plpl,pl;ml,mimi) = C (^^ )P ?;m 2 2) m|,m?) , (Bll) 
C (pl,pl,pl-,ml,ml,n4) = C (pj, p\, p\\ m 2 , m 2 , m 2 ) . (B12) 

(2) Do-function 

n/2222 2 2 2 2\ 

A) (Pl,P 2 ,P 3 > P45 s 12, S23; Win m 2 ? m 3> m i) 

= D (pl } pl J pl } p 2 1 -,s 2 3,s 12 -,ml,ml,ml,m 2 1 ) , (B13) 

n/2222 222 2\ 

A) (Pi,P 2 )P3)P4; si2,s 23 ; m 1 ,m 2 ,m 3 ,m 4 J 

= A) (piP3^2^i; S 12,S23;^?,m4>m3> m 2) • ( B14 ) 

APPENDIX C: FORM FACTORS OF THE VIRTUAL CORRECTIONS 

• The form factors of the triangle loop V\ 

1 

h 



f 2 Vl = ^{2C Q mth + 2m t C u ti \, (CI) 



/ 6 Vl = U 80 ^ + 8m 2 C n - 4m 2 C 21 + 8C 24 |, (C2) 
= ^ {2Cb (2m 2 + + 2 (4m 2 + C u + 2t x C l2 

+2m 2 t C 21 + 2t x C 2Z + 4C 24 J , (C3) 

fll = i|4m t Cn+4m t C 21 |, (C4) 

where the arguments of the scalar function and tensor coefficients are 
{{-Ptf , P% (Pt - P 9 f ; 0, m 2 , m 2 ) = (m 2 , 0, t; 0, m 2 , m\). 

The form factors of the triangle loop V 2 

ft* = LUcomth + BmtdihX, (C5) 



fe 2 = ^|l2C m 2 + 12C 21 m 2 + 2 (12m 2 + ^) C n + 4t 1 Ci 2 + 8tiC 23 + 24C 24 

+e (-8C m 2 - 16Cnm 2 - 8C 2 im 2 - At x C 12 - 4t 1 C 23 - 16C 24 ) J, (C6) 
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f% = ^6C m 2 + 2C 21 m 2 + 2(4m 2 + t 1 )C 11 + t 1 C 12 

+2hC 23 + 12C 24 - 8eC 2 4 \ , (C7) 



fii = t~ <! -4m t C n - 4m t C 21 + e (4C m t + 8Cnm t + 4C 21 m t ) } , (C8) 



where the arguments of the scalar function and tensor coefficients are 
{(-P2) 2 , V% (Pt ~ P g f ; ml 0, 0) = K, 0, t; m 2 , 0, 0). 

The form factors of the triangle loop V3: 

= ^{4C (m\ - m 2 w + t ± ) + 4 (m 2 - m 2 w + t ± ) C n + 4 (m 2 + t x ) C 12 



+4m 2 w C 22 + A(m 2 t -m 2 w + t l ) C 23 + 8C 24 - 16C 24 e| ; 



(C9) 



fg 3 = ^|-8m t C 12 + 8em t C 23 -8m t C 23 |, (C10) 
= ^|2C (m 2 - m 2 ^ + t ± ) + 2 (m 2 - m 2 w + C n + 2 (m 2 + ^) C 12 
+2m 2 w C 22 + 2 (m 2 - m 2 ^ + C 23 + 4C 24 - 8C 24 e j, (Cll) 
fii = U-4m t C 12 -4m t C 23 X, (C12) 



where the arguments of the scalar function and tensor coefficients are 
{pI {-Pw? , {Pw ~ Pbf ; 0, 0, m\) = (0, m 2 w , t; 0, 0, m 2 ). 

• The form factors of the triangle loop V{ ( fj 1 = F^ 1 / s) 

fl ( = (8C 24 - 16C 24 e)/s, (C13) 
= 4C + 4C n + 4C 12 + 4C 23 + 8C 24 /s + (-4C 12 - AC 23 - 16C 24 /s) e, (C14) 
/ff = -4C - 4C n - 4C 12 - 4C 23 - 8C 24 /s + (4C 12 + 4C 23 + 16C 24 /s) e, (C15) 
f x i = 2C + 2C n + 2C 12 + 2C 23 + 4C 24 /s + (-2C 12 - 2C 23 - 8C 24 /s) e, (C16) 

where the arguments of the scalar function and tensor coefficients are 



^ 2 ,p 2 ,(^ + P 9 ) 2 ,0,0,0) = (0,0,s; 0,0,0). 
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The form factors of the triangle loop V 2 . 

f¥ = 2C U + AC 12 + 8C 23 + 2AC 24 /s + e (-4C 12 - 4C 23 - 16C 24 /s) , (C17) 

ff = AC U + 2C 12 + 4C 23 + 24C 24 /s - 16eC 24 /s, (C18) 

/n = -4C n - 2C 12 - 4C 23 - 24C 24 /s + 16eC 24 /s, (C19) 

f\i = 2C n + C12 + 2C 23 + 12C 24 /s - 8eC 24 /s, (C20) 

where the arguments of the scalar function and tensor coefficients are 
(pl P% (Pb + Pgf , o, 0, 0) = (0, 0, s; 0, 0, 0). 

The form factors of the triangle loop V 3 L . 

fi* = -^AComt + Amt^, (C21) 
ff = -2C m t - 2m t C 11 , (C22) 
/T 3 ' = ~{ 4C7 o W - + s) + 4 (2m 2 - + s) C u + (4s - Am 2 ) C 12 

+Am 2 t C 21 + Am 2 w C 22 -A(m 2 t +m 2 w - s) C 23 + 8C 24 - 16C 24 e| , (C23) 

= - s ^AC (m 2 t - m 2 w ) + (8m 2 t - Am 2 w ) C u - Am 2 t C 12 

+Am 2 t C 21 + Am 2 w C 22 -A(m 2 t + m 2 w ) C 23 + 8C 24 - 16C 24 e| , (C24) 



= 1 < -BmCu + 8m t C 12 - 8m t C 21 + 8m t C 23 }, (C25) 



1 



s 



fS = - S {- 4C ° K -m 2 w + s)-A (2m 2 t - m 2 w + s) C n + 4 (m 2 t - s) C l2 

-Am 2 C 21 - Am 2 w C 22 + A (m 2 t + m 2 w - s) C 23 - 8C 24 + 16C 24 e|, (C26) 
= ^2C (m 2 t -m 2 w + s)+2 (2m 2 - m 2 w + s) C u + (2s - 2m 2 ) C 12 

+2m 2 t C 21 + 2m^C 22 - 2 (m 2 + m 2 ^ - s) C 23 + 4C 24 - 8C 24 e| , (C27) 



fii = ^ I -4m t C n + Am t C 12 - Am t C 21 + Am t C 23 j , (C28) 

where the arguments of the scalar function and tensor coefficients are 
{{-Ptf , {-Pwf , (Pt + Pwf , 0, m 2 , 0) = (m 2 , m 2 w , s; 0, m 2 , 0). 



1 
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The form factors of the box loop B\\ 



ft 



fB 1 
J2 



J4 



rBi 
J5 



Je 



m t |4 (2m 2 t + h) D 21 + 4 (m 2 + m 2 w + h) D 23 + 4t 1 L> 24 

-4 (3m 2 + m 2 w + 2t,) D 25 - U X D 2% + 16D 27 + AD 31 m 2 - Am 2 w D 33 

+4t 1 L> 34 - 4 (2m 2 + mjv - s) L> 35 + 4 (m 2 + 2m 2 w - s) D 37 

+4 (s + ix) D39 - 4 (s + 2h) D 3W + I6D311 - I6D313 1, 



77it|-2s£)ii + 2sL>i 3 - 2 (2m 2 + D 21 - 2 (m 2 + m^ + fx) L> 23 
-2txD 24 + 2 (3m 2 + m 2 w + 2ti) L> 25 + 2txD 26 - 8^27 - 2D 31 m 2 



+2m 2 w D 33 - 2t lJ D 34 + 2 (2m 2 + 



>) D 35 - 2 (m 2 + 2m 2 w - s) D 37 



-2 (s + h) D 39 + 2 (s + 2t ± ) D 3W - 12D 311 + 12D 313 



+e 



-2D 31 m 2 - 2 (2m 2 + h) D 21 - 2 (m 2 + m 2 w + h) D 23 - 2t 1 L> 24 
+2 (3m 2 + m 2 w + 2t x ) D 25 + 2t lJ D 26 - 4D 27 + 2m 2 w D 33 - 2t 1 D 34 
+2 (2m 2 + m 2 w - s) D 35 m t - 2 (m 2 + 2m 2 v - s) D 37 m t - 2 (s + D 39 m t 



+2 (s + 2ti) D 3 io™t - 8£> 3 nmt + 8D 3 i 3 "it 



+ e 2 (4£> 3 ii - 4L>. 



313, 



-8m 2 L>ii - ttiD 12 + 4 (2m 2 + h) D 13 - Am 2 D 21 - 4ti_D 24 + 4m 2 L> 25 



+4tiD 26 - 8L> 27 - 8D 3 i 3 + e 



-4,035m 2 + 4ti^ 23 - 4tiD 25 + 8D 27 



-Am 2 w D 33 + 4 (m 2 + m^ - s) D 37 + 4 (s + h) D 39 - At x D 3W 



-16D 



313 



+ 8D 3 i 3 e 2 , 



= -4*i£>i 2 + 4<iDi 3 - 4m 2 J D 21 - Am 2 w D 23 - 4hD 24 + 4 (m 2 + m^) £> 



25 



+4^1)26 - 8D 27 - 8D311 + e 



(C29) 



(C30) 



(C31) 



-8D 21 m 2 - AD 31 m 2 + 4t 1 £> 22 + (8m 2 - 4^) £> 24 
+4 (m 2 + m^) D 25 - 4 (m 2 + m 2 w ) D 26 + 8D 27 + 4 (m 2 - L> 34 
+4 (m 2 + m^ - s) D 35 + 4t 1 D 36 - Am 2 w D 37 - 4 (s + *x) L> 38 + Am 2 w D 39 
-4 (m 2 + m^ - 2s - ^) D 310 - 1QD 3U + I6D312 

+8D3i2(8 J D 3 i 1 -8Am)e 2 ( C32 ) 
AD (m 2 - m 2 w + s + + (12m 2 + 8 (-m 2 ^ + s + h)) D u + (4*i - 4s) D 12 

+ (-8m 2 + Am 2 w - 8t ± ) D 13 + 4 (3m 2 - m 2 ^ + s + D 2 i + 4 (m 2 + m 2 ^ + D 23 

+ (8*i - 4s) D 24 - 8 (2m 2 + D 25 - 8^L> 26 + 16D 27 + 4m 2 D 31 - 4m 2 4 ,D 33 
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+4t 1 £>34 - 4 (2m 2 t +m 2 w -s) D 35 + 4 [m\ + 2m 2 w - s) D 37 + 4 (s + t ± ) D 39 
-4 + 2h) D 3l0 + I6D311 - I6D313 

+ (AD 21 m 2 + AD 23 m 2 t - 4 (2m 2 + s) D 25 + AsD 26 - 8D 27 ) e + 8D 27 e 2 (C33) 
/ 7 Bl = 8m t jo^ - 2/^5 + D 26 + D 39 - D 3W + e (D 3W - D 39 ) j (C34) 

fi 1 = Sm t ^D 25 - D 2% - D 35 + D 37 - D 39 + D 310 + e (D 35 - D 37 + D39 - Ano) |, (C35) 

/ 9 Bl = 8m J -D u + 0> 12 - 2D 21 - D 23 + 3^25 - D u - D 39 + 2D 



.; - OY/i t ^— u\\ -1- 1^12 — ^-^21 — ^23 T <3-t^25 — -^34 ~ -^39 ~r ^-1^310 

+e (^24 - ,026 + D 34 + D 39 - 2D 310 ) ) , (C36) 



/fo 1 = 8m t j-£> 12 + 0> 13 - 0> 3 i + ^34 + 2D 35 - D 37 + D 39 - 2D 310 

+e (L>2i - D24 - D 25 + D26 + An - D u - 2D 35 + D 37 - D 39 + 2 A10) J , (C37) 
/n 1 = 40)21^ - 80) 2 4m 2 - 4t!0)22 - 4 (m 2 + 0> 23 + 4 (m 2 - m 2 w + s + 0> 25 
+4 (m t 2 + - s + h) D 26 - lQeD 27 + Am 2 w D 33 - AD 34 m 2 + 4D 35 m 2 
-AhDsz - 4 (m 2 + m 2 w -s) D 37 + 4 (s + t x ) D 38 - 4 (m 2 w + s + D39 
+4 (m 2 + m^ - s + h) D 310 - 1QD 312 + I6D313 (C38) 
f$ = 2D (m 2 - m 2 w + s + + 2 (2m 2 - m^ + s + D n - 2m 2 0> 13 - 2m 2 0> 2 i 
+2t 1 D 22 + 2 (m 2 + fx) 0) 23 + 4m 2 0>24 - 2 (m 2 - m 2 w + s + t x ) D 25 
-2 (m 2 + m 2 w -s + h) D 26 + AD 27 - 2m 2 w D 33 + 2m 2 D 34 - 2m 2 D 35 
+2^36 + 2 (m 2 + m 2 w - s) D 37 - 2 (s + ^) As + 2 (m^ + s + h) D 39 
-2 (m 2 + m 2 w -s + t ± ) D 310 + 12D 312 - 12D 313 + (AD 27 - AD 312 + 4D 3 i 3 ) e 2 
+e -2D 21 m 2 + AD 24 m 2 + 2D u m 2 - 2D 35 m 2 + 2t x D 22 + 2 (m 2 + ^) 0> 23 
-2 (m 2 -m 2 w + s + h) D 25 - 2 (m 2 + m 2 w - s + D 26 + 8D 27 - 2m^D 33 
+2t!0)3 6 + 2 (m 2 + m^ - s) 0> 37 - 2 (s + h) D 38 + 2 (m 2 w + s + h) D 39 
-2 (m 2 + m 2 w -s + t ± ) D 310 + 8D 312 - 8D 313 (C39) 
fw 1 = -4m t D n +Am t D 12 + 4m t D 24 -4m t D 26 + e(4m t D 24 -4m t D 2(i ) (C40) 
f$ = 4m t D u - Am t D l3 - Am t D 25 + e (4m t D 23 - Am t D 25 ) + 4m i D 26 (C41) 
fEi = -Am t D l2 + 4m t 0' 1 3 + Am t D 21 - Am t D 2A - Am t D 25 + Am t D 26 

+e (4m t O) 2 i - 4m t O> 24 - 4m t O) 25 + 4m t O> 26 ) (C42) 

34 



J 16 



Jl7 



il8 ~~ 



J 19 



i20 



Am t D 12 - Am t Di 3 - Am t D 2 i + Am t D 24 + ^m t D 25 

+e (-Am t D 21 - Am t D 23 + 8m t D 25 ) - Am t D 26 (C43) 
~8D 23 + 8^26 + 8D38 - 8D39 + e (8^23 - %D 2& - 8D 38 + 8D39) (C44) 
-8D25 + 8D 26 - 8D 37 - 8D 38 + 8D 39 

+e (-8D 25 + 8D26 + 8D 37 + 8D38 - 8D 3g - 8D 3l0 ) + 8D310 (C45) 
16D 12 - 16D 13 + 8D22 + 8D23 + I6D24 - 8D25 - 24^26 + 8D36 - 8D 38 + 8D 39 
-8D 3W + e (-8D22 - 8D23 + 8D25 + 8D 26 - 8D 36 + 8D 3 8 - 8D 39 + 8D310) (C46) 
-8D22 + 8D 24 - 8D 25 + 8D26 + 8D34 - 8D35 - 8D36 + 8D 37 + 8D 38 - 8D 39 



+e 



8D 



22 



8D 24 + 8D25 - 8D 26 - 8D34 + 8D35 + 8D. 



36 



-8D37 - 8D38 + SD. 



39 



(C47) 



where the arguments of the scalar function and tensor coefficients are 
{(-Ptf , P 2 g , {-Pwf , P\\ t, u; 0, m 2 , m 2 , 0) = (m 2 , 0, m 2 w , 0; t, u; 0, m 2 , m 2 , 0). 

• The form factors of the box loop B 2 

= -AD 21 m 3 t - 4 (m 2 - m 2 w + ij) D xl m t + 4 (m 2 - m 2 w - s + D 12 m t + 4sD 13 m t 
-Am 2 w D 22 m t + 4 (m 2 + m 2 ^ - s) D 24 m t - At x D 2 ^m t + A(s + t 1 ) D 26 m t - 8D 27 m t 
-8D 311 m t + 8D 312 m t + e 2 (8m t D 3U - 8m t D 312 ) + e -AD 31 m? t 
-4 (m 2 - m 2 ^ + s) D 21 m t + 4 (-m 2 + m 2 w + s) -D 2 2«^ + 8 (m 2 - m 2 ^) -D 2 4^i 
+4sD 25 m t - AsD 26 m t + 8D 27 m t + Am 2 w D 32 m t + 4 (2m 2 + m 2 w - s) D 3A m t 
-AtiD^mt - 4 (m 2 + 2m 2 w - s) D 36 m t - 4 (s + t x ) D 38 m t + 4 (s + 2ti) £310^ 



-16D 311 m t + 16D 3 i 2 m t 



(C48) 



/ 2 B2 = 2D 31 m\ + 2 (m 2 - m 2 w + h) D X \m t + 2 (-m 2 + m 2 ^ + s - h) D 12 m t - 2sD 13 m t 



+2 (2m 2 - m 2 ^ + s) D 21 m t + 2 (m 2 - s) D 22 m t + 2 (-3m 2 + m 2 ^ + s) D 2A m t 



+2 (£i - s) -D 2 5^t - 2t 1 D 26 m t + 4D 27 m t - 2m 2 w D 32 m t - 2 (2m 2 + 



m 



w 



i) D u m t 



+2t 1J D 35 m t + 2 (m 2 + 2m^ - s) D m m t + 2 (s + t x ) D 38 m t - 2 (s + 2t x ) Dgio^t 



+12D 311 m t - \2D 312 m t + e 



2^31^? + 2 ( 



m 



w + s ) D 2im 



+2 (m 2 - m 2 ^ - s) -D 2 2^t - 2sD 25 m t + 2sD 26 m t - 4D 27 m t - 2m 2 w D 32 m t 
-2 (2m 2 + m 2 w - s) D 3A m t + 2tiD 35 m t + 2 (m 2 + 2m 2 w - s) D 36 m t 
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+2 (s + h) D 38 m t - 2 (s + 2ti) £310^ + 8^311^ - 8D 312 m t 



+ [Am t m 2 w - Am\ ) -°24 



+ e 2 (Am t D 312 - 4m t D 311 ) , 



(C49) 



AD ( 



m. 



m 



w 



+ s + ti) + 4 (2m 



m 



+ s + ti) D u + (8s - Am 2 w ) D 12 



+ {Ah - As) D13 + 4m 2 D 2 i + 4 (-m 2 + + s) D 22 + 4 (m 2 — 2m^ + 2s) L> 24 
+4t lJ D 25 - 8sD 26 + 1QD 27 + 4m^£> 32 + 4m 2 J D 34 - 4 (m 2 + 772^-5) D 36 
-4 (s + ix) D 3 8 + A^Dsw + I6D312 

+ (-AD 21 m 2 t - AD 22 m 2 t + 8D 24 m 2 - AtiD 25 + At x D 2% - 8D 27 ) e + 8D 27 e 2 , (C50) 
AD (m 2 - m 2 w + s + + (12m 2 + 8 (-m 2 ^ + s + fi)) D u 



+4 (-2m 2 + m 2 w + s - h) D l2 - 8sL>i 3 + 4 (2m 2 - m 2 ^ + s + ti) D 2 i + 4m 2 v D 22 
+4sL> 23 - 4 (2m 2 - 2s + t±) D 24 - A (m 2 - m 2 w + 2s) D 25 + A (m 2 - m 2 w - 2s) L> 26 
+16D 27 + 4m 2 D 3 i - 4 (m 2 + m 2 w - s) D 3 4 + (4*i - 4m 2 ) D35 + 4m 2 ^ J D 3 6 - 4tiD 37 
-Am 2 w D 38 + 4 (s + h) D39 + 4 (m 2 + m^ - 2s - D 310 + I6D311 - 16£> 3 i3 
+ (-AD 21 m 2 - Am 2 w D 22 + A (m 2 + m 2 w ) D 24 - AhD 25 + At\D 26 - 8D 27 ) e 
+8D 27 e 2 , (C51) 
-AsD 12 + 4sD 13 - AsD 2A + 4sD 26 - 8D 27 - 8D311 + 8D312 + (8D311 - 8D312) e 2 
+e -AD 31 m 2 + 4s£> 22 - 4sL> 24 + 8L> 27 + Am 2 w D 32 + A (2m 2 + m 2 w - s) D 34 

-4ti£>36 - 4 (m 2 + 2m 2 w - s) D 36 - 4 (s + h) D 38 + 4 (s + 2ti) D310 

-16^311 + 16^312 , 

-8m t D l2 + 8m t D 13 - 8m t D 24 - 8m t D 25 + 16m t L> 26 + 8m t D 38 - 8m t D 3W 
+e (8m t D 22 - 8m t D 24 + 8m t D 25 - 8m t D 26 - 8m t D 38 + 8m t D 3l0 ) , 
-8m t Dii + 16m t D 12 - 8m t D 13 - 8m t D 21 + 8m t D 22 + 8m t D 24: + 8m t D 25 



(C52) 



(C53) 



-16m t £) 26 - 8m t D 3A + 8m t D 36 - 8m t D 38 + 8m t D 310 + e 



8m t D 21 - 8m t D 22 



-8m t D 25 + 8m t D 26 + 8m t D 3A - 8m t D 36 + 8m t D 38 - 8m t D 310 , 
-8m t D 2l + 8m t D 24 + 8m t D 25 - 8m t D 26 - 8m t D 35 - 8m t D 38 + 16m t D 3 i 
+e (-8m t D 21 - 8m t D 22 + 16m t D 24 + 8m t D 35 + 8m t D 38 - 16m t D 3W ) , 
-8m t D 22 + 8m t D 24 - 8m t D 25 + 8m t D 26 - 8m t D 31 + 16m t D 34 + 8m t D 35 



(C54) 



(C55) 



-8m t D 36 + 8m t D 38 - 16m t D 310 + e 



8m t D 2 i + 8m t D 22 - \<6m t D 2i 
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+8m t D 31 - 16m t D u - 8m t D 35 + 8m t D 36 - 8m t D 38 + 16m t D 310 
(-8D 27 + 8D 312 - 8D 313 ) e 2 

+e 4D 21 m 2 - 8D 2A m\ - 4D 34 m 2 + 4:D 35 m 2 + 4 (m\ - s) D 22 - 4sD 23 
+4 (m 2 -mly + s + h) D 25 - 4 (m 2 - m 2 w - s + h) D 26 - Am 2 w D 32 
+4 (m 2 + m 2 w -s) D 36 + tt x D 31 + 4 (m^ + s + h) D 38 - 4 (s + t x ) 7J 39 
-4 (m t 2 + - s + h) D 3W - 1QD 312 + I6D313 
-AD (m 2 -m 2 w + s + h) - 4 (2m 2 - m 2 w + s + h) D n + Am 2 D 12 
—8D 27 — 8D312 + 8-D313 

-2D 21 m 2 + 4 J D 24 m 2 + 2D u m 2 - 2D 35 m 2 + (2s - 2m 2 ) D 22 + 2sD 23 



(C56) 



(C57) 



-2( 



m 



w 



+ s + tx) D 25 + 2(m 2 -m 2 w -s + h) D 2% + 8D 27 + 2m 2 w D : 



32 



-2 (m 2 + m<V - s) D 3% - 2t x D 37 - 2 (m 2 w + s + *J 7J 3 8 + 2 (s + ^) £>. 



39 



+2 (m 2 + - s + tx) £> 310 + 8D 312 - 8D 



313 



+2D (m 2 - m 2 w + s + h) + 2 (2m 2 -m^ + s + h) D n - 2m 2 D 12 - 2m 2 D 21 



+ (2s - 2m 2 ) D 22 + 2sD 23 + Am 2 D 24 - 2 ( 
+2 (m 2 - m 2 w - 



m 



w 



+ s + h) D 



25 



fi) D 26 + 4D 27 + 2m 2 w D 32 + 2m 2 D u - 2m 2 D 35 



-2 (m 2 + m 2 w -s) D 36 - 2t t D 37 - 2 (m 2 w + s + t ± ) D 38 + 2 (s + t ± ) D 39 
+2 (m 2 + m 2 w -s + t x ) D 3W + 12D 312 - 12D 313 
+ (AD 27 - AD 312 + AD 313 ) e 2 

-4m t D u + Am t Di 2 + 4m t D 25 - Am t D 2% + e (4m t D 25 - Am t D 2& ) 
Am t Du - 8m t D 12 + 4m t D 13 - 4m t D 2A + e (4m t D 22 - 4m t D 24 ) + 4m t D 26 
4m t D 21 - 4m t D 24 - 4m t D 25 + 4m t D 26 
+e (4m t D 21 - Am t D 2A - 4m t D 25 + Am t D 2% ) 
-Am t D 21 + 4m t D 24 + e (-4m t D 21 - Am t D 22 + 8m t D 24 ) + 4m t D 25 - Am t D 2% (C62) 
8D 23 - 8D26 - 8D38 + e (-8D23 + 8D 26 + 8D 38 - 8D39) + 8D39 (C63) 

-I6D12 + I6D13 - 8^22 " 8^23 " 16^24 + 8^25 + 24D 26 " 8^36 + 8£ 3 8 - 8^39 

+e (8D 22 + 8D23 - 8D25 - 81)26 + 8D36 - 8D38 + 8D39 - 8D310) + 8D310 (C64) 
8D 25 - 8D 26 + 8D 37 + 8D 38 - 8D 39 - 8D 310 



(C58) 
(C59) 
(C60) 

(C61) 
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+e (8D 25 - 8D 26 - 8D 37 - 8D 38 + 8D 39 + 8D 310 ) (C65) 
f^ = 8D 22 - 8D 2A + 81)25 - 8D 26 - 8D 34 + 8D 35 + 8D 36 - 8D 37 - 8D 38 + 8D 39 



+e 



-8£> 22 + 8D 24 - 8D 25 + 8D 26 + 8D 34 - 8D 35 - 8D 36 



+8D 37 + 8^38 - 8D 39 



(C66) 



where the arguments of the scalar function and tensor coefficients are 
{{-Ptf , {-Pwf , P 2 g , vl\ s, u; 0, m 2 , 0, 0) = (m 2 , m 2 w , 0, 0; s, u; 0, m 2 , 0, 0). 

The form factors of the box loop B 3 : 

fj* = -AD 21 m 3 t - 2D 31 m 3 t + At x D xx m t - 4^D 12 m t - 2 (m 2 t + m 2 w - s) D 22 m t 
+2 (3m 2 + m 2 w - s) D 24 m t + 2 (m 2 - m 2 ^ + s + t 4 ) D 25 m t 
-2 (m 2 - m 2 w + s + ti) D 26 m t - \2D 27 m t + 2m 2 v D 32 m t 
+2 (2m 2 + m 2 w - s) D u m t + 2 (m 2 - m 2 w + s + t x ) D 35 m t 
-2 (m 2 + 2m^ - s) D 36 m t + 2 (m 2 - m 2 w + ^) D 38 m t 
-2 (2m 2 - 2m^ + s + 2ti) Ajio™* - 4D 311 m t + 4 J D 312 m t 

+e(4m tJ D 3 i2-4m tJ D 3 ii) (C67) 
/.f 3 = 2m t D 22 m 2 w + m t (s- 2ti) D u -m t (s- 2t ± ) D 12 + m t (m 2 + m 2 ^ - s) D 21 
+m t (-m 2 - 3m 2 w + s) L> 2 4 - m t (m 2 t - m 2 w + ti) D 25 

+m t (m 2 t - m 2 w + t x ) _D 2 6 + Qm t D 27 (C68) 
/f 3 = -ADot, +A(m 2 - t,) D n - 2 (3m 2 - m 2 w + s) D 12 + 2 (m 2 - m 2 w + s + h) D l3 
+2m 2 D 21 + 2 (m 2 - s) D 22 - 2 (m 2 - m 2 w + s + ti) D 23 - 2 (3m 2 - m 2 w + s) D 24 
+2 (m 2 -m 2 w + s- t x ) D 25 + 2 (m 2 - m 2 w + 2 (s + ^)) Dae - 2m 2 yJ D 32 
-2m 2 D 34 + 2m 2 D 35 + 2 (m 2 + m^ - s) D 36 - 2 (m 2 - m 2 w + s + fx) D 37 
-2 (m 2 - 2m 2 w + h) D 38 + 2 (m 2 - m^ + ij L> 39 + (-4m^ + 4s + 2ti) .Dsio 
-4.D312 + 4D 3 i3 + e AD 3A m 2 t - AD 35 m 2 + 4sD 2 2 + 4sD 23 - 8sD 26 
-4D 27 + 4m 2 w D 32 - 4 (m 2 + m 2 w - s) D 36 + 4 (m 2 - m 2 w + s + ^) D 37 
+4 (m 2 - 2m^ + h) D 38 - 4 (m 2 - m^ + t x ) D39 + (8m^ - 4 (2s + h)) D 310 
+20D 312 -20D 313 + (8D 313 - 8D 312 ) e 2 (C69) 
/f 3 = 2 J D 31 m t 2 - 4 J D t 1 + 4 (m 2 - 2^) - 2 (m 2 + m^ - s - 2*i) D 12 
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-2 (m 2 - m 2 w + s + h) D 13 + (lOm 2 - 4^) D 21 + 6m 2 w D 22 

+4 {-2m 2 - 2m 2 w + s + t ± ) D 24 - 4 (2m 2 - 2m 2 w + s + 2t ± ) D 25 

+6 (m 2 -m 2 w + ti) D 26 + 28D 27 - 2 (m 2 + m 2 w - s) D 34 

-2 (m 2 t - m 2 w + s + ti) D 35 + 2m 2 v D 36 + 2 (m\ - m 2 w + ti) D 310 

+e {-20D 27 - 12D 311 ) + 2OA311 (C70) 

ft = -2D s - 2 (m 2 -m 2 w + 2s) D n + 2 (m 2 - m 2 w + 2s) D 12 + 2sD 13 
-2 (2m 2 -m 2 w + s) D 21 + (2s - 2m 2 ) D 22 + 2 (3m 2 - m 2 w + s) D M 
+2 (m 2 -m 2 w + U) D 25 - 2 (m 2 - m 2 w + s + h) D 26 - 2m 2 D 31 + 2m 2 w D 32 
+2 (2m 2 + m 2 w -s) D 34 + 2(m 2 -m 2 w + s + h) D 35 - 2 (m 2 + 2m 2 w - s) D 36 
+2 (m\ -m 2 w + h) D 38 - 2 (2m 2 - 2m 2 w + s + 2ti) D 310 - AD 3ll + AD 312 
+e AD 31 m 2 - AsD 22 + AsD 2A - AsD 25 + AsD 26 - AD 27 - Am 2 w D 32 
-A (2m 2 + m 2 w - s) D 34 - 4 (m 2 - m 2 w + s + h) D 35 
+4 (m\ + 2m 2 w - s) D 36 - 4 (m\ - m 2 w + h) D 38 
+4 (2m 2 - 2m 2 w + s + 2ti) A310 + 20D 311 - 20D 312 
+ (8D 312 -8D 311 )e 2 (C71) 

ft = 8m t D 12 - 8m t D 13 + 8m t D M - 8m t D 25 + 8m t D 37 + 8m t D 38 
-8m t D 3g - 8m t D 310 + e I -8m t D 22 + 8m t D 24 - 8m t D 25 



+8m t D 26 - 8m t D 37 - 8m t D 38 + 8m t D 39 + 8m t D 310 j (C72) 
ft = -8m t D 12 + 8m t D 13 - 8m t D 22 + 8m t D 26 + 8m t D 34 - 8m t D 35 

-8m t D 36 + 8m t D 3W + e I 8m t D 22 - 8m t D 2A + 8m t D 25 - 8m t D 26 



-8m t D 34 + 8m t D 35 + 8m t D 36 - 8m t D 3W j (C73) 
ft = 8m t D 21 - 8m t D 2A - 8m t D 35 - 8m t D 38 + 16m t D 3W 

+e (8m t D 21 + 8m t D 22 - lQm t D 24 + 8m t D 35 + 8m t D 38 - lQm t D 310 ) (C74) 
ft = 8m t D 22 - 8m t D 2A + 8m t D 31 - 16m t D 34 + 8m t D 3% 

+e (-8m t D 21 - 8m t D 22 + l6m t D 24 - 8m t D 31 + 16m t D u - 8m t D 36 ) (C75) 
ft = -6D 21 m 2 - 2D u m 2 + ADoh + (At Y - Am 2 ) D u + 2 (m\ + m 2 w - 3s) D 12 

+2 (m 2 -m 2 w + 3s + t ± ) D 13 + 2 (m 2 - 2m 2 w - s) D 22 + 2 (m 2 + 3m 2 w - 3s) D 24 
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+6 (m 2 -rn^ + s + h) D 25 + {~Am 2 t + Am 2 w + 2s - 4^) D 2% - 8D 27 - 2m 2 w D 32 
+2 (m 2 + m 2 w - s) D m - 2 [m\ - m 2 w + ti) D 38 + 2 [m\ - m 2 w + s + h) D 310 

+e(12D 27 -AD 312 )-AD 312 (C76) 
f$ = 5D 21 m 2 - 2£> *i + 2 (m 2 - D u + (-m 2 - m 2 w + 4s) 0> 12 

+ (-m? + - 3s - ii) £) 13 + 3m^O> 2 2 - 4 (m 2 + m 2 w - s) D 24 

-5 (m 2 -m 2 w + s + t x ) D 25 + A (m 2 - m 2 w + ti) D 26 - UeD 27 + QD 27 (C77) 

ffi = Am t D n - 4m t D 12 + 8m t D 25 - 8m t D 2% (C78) 
= -4m t D u + 8m t D l2 - 4m t D 13 + 4m t D 2A - 4m t D 25 

+e (-Am t D 22 + Am t D 24 - Am t D 25 + Am t D 26 ) (C79) 

Af = -4m t D n + 4m t D 12 - 8m tJ D 21 + 8m t D 24 (C80) 

At = Am t D 21 - Am t D 2A + e (Am t D 21 + Am t D 22 - 8m t D 24 ) (C81) 

fi 7 3 = 8D 23 - 8D 26 - 8D 38 + e (-81^3 + 8D26 + 8^33 - 8D39) + 8D 39 (C82) 
= I6A2 - 16^13 + 8D22 + 16^24 - 1QD 25 - 8D 2% + 8D36 - 8D310 

+e (-8D22 + 8D26 - 8/^36 + 8D310) (C83) 

At = 8D25 - 8D26 + 8D38 - 8D310 + e (8D25 - 8D26 - 8D38 + 8D310) (C84) 

At = -8,022 + 8,024 + 8,034-8,036 + 6(8.022-8,024-8,034 + 8,036), (C85) 

where the arguments of the scalar function and tensor coefficients are 
{{-Ptf , {-Pwf , Pi P 2 g \ s, t; 0, m 2 , 0, 0) = (mf , m 2 w , 0, 0; s, t; 0, m 2 , 0, 0). 

Bubble correction Si. 

= ^4B m t + 2B 1 m t + e(-2B m t -2B 1 m t )^, (C86) 

fe 3 = ^l6B m 2 + 4B 1 (2m 2 + t 1 )+e(-8B m 2 -4B 1 (2m 2 + t 1 ))^, (C87) 

Af = ^8B m 2 t + 2B X (2m t 2 + t,) + e (-45 m? - 25 x (2m? + ^)) j, (C88) 

where the arguments of the scalar function and tensor coefficients are 
{(Pg -Pt) 2 ,m 2 ,0) = (t,m 2 ,0). 
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Bubble correction S2 : 
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2B ± e} , 



(C89) 
(C90) 
(C91) 
(C92) 

where the arguments of the scalar function and tensor coefficients are 
(( Pg + p b ) 2 ,0,0) = (s,0,0). 

APPENDIX D: DIVERGENCES OF THE FORM FACTORS 



In this section we list out the divergent pieces of the form factors and and we further 
distinguish the UV divergence and the IR divergence. The former is written inside the 
square brackets, i.e. [•••]. 

(1) The triangle loop V\ only gives rise to the UV divergences, 



Vi 
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V12 = T 
ti 
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(2) The triangle loop V 2 exhibits both the UV and IR divergences, 
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(3) The triangle loop V 3 gives rise to both the UV and IR divergences, 
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(4) The triangle loop V[ give rise to both the UV and IR divergences, 
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(5) The triangle loop V 2 ' gives rise to both the UV and IR divergences, 
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(6) The triangle loop V% contains only the UV divergences: 
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(7) The box loops #1,2,3 only contain the IR divergences. B x gives rise 
divergences: 
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(8) B 2 gives rise to the following divergences: 
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(9) B 3 gives rise to the following divergences: 
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(10) The bubble loops S 1:2 only contain the UV divergences. Si gives rise to the 
divergences: 
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